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relativistic covariance. The sum rules can be used for isospin-asymmetric nu- 



the isospin density to weighted sums over states with fixed velocity relative 
to the nuclear matter ground state. 
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Sum rules have been an important tool in theoretical physics for a very long time. Their 
importance stems in large measure from the fact that they allow one to make concrete 
predictions concerning sums of matrix elements even when one does not have, or cannot solve 
the theory for the dynamics of the system. One interesting and important system for which 
the underlying dynamics is known but cannot presently be solved is the nuclear medium. 
While we are quite onfident that the underlying dynamics for the strong interactions is QCD, 
it is quite unlikely that it will prove to be a tractable task to calculate nuclear properties from 
QCD any time in the forseeable future. This raises an obvious question, namely whether one 
can use known properties of QCD, such as current algebra, to derive sum rules for nuclear 
properties. 

In this work we will not discuss sum rules for finite nuclei. Instead, we concentrate on 
infinite nuclear matter with arbitrary baryon density and arbitrary isospin density. We do 
this for several reasons. One is practical: if the system is translationally invariant (and 
hence necessarily infinite) momentum is a good quantum number and can be used to label 
excitations. A second reason is that interiors of large nuclei are well approximated by nuclear 
matter and learning about nuclear matter gives a more general insight into nuclear systems 
than the study of individual nuclei. A third reason is that the cores of neutron stars are 
essentially isospin-asymmetric nuclear matter (neutron matter). Finally, infinite nuclear 
matter is theoretically interesting in its own right. There is obviously one major downside 
to studying infinite nuclear matter: the connection between sum rules for infinite nuclear 
systems and experimental observables for real nuclei is not immediately clear. 

In previous works |]l],0] we have studied the chiral properties of nuclear matter. That 
is, we examined the behavior of modes as the quark mass is formally taken to zero. We 
showed that certain zero-momentum modes had to go to zero energy as the (current) quark 
masses were taken to zero. In the process of deriving this result we used current algebra and 
completeness to derive sum rules for the squares of the matrix elements of the divergence of 
the axial current between zero momentum states. Analogous sum rules are also implicit in 
the work of Lutz, Steiner and Weise P|. 



Now, in Refs. [l|,^ the central stress was on the fact that these sum rules were necessarily 
saturated by modes whose energy went to zero as the quark masses did. However, it is worth 
observing that these sum rules in no way depend on the chiral limit. As sum rules, they 
hold for any value of the quark mass. This is significant for two reasons. The first is that 
compared to typical nuclear mass scales, the effects of finite up and down quark masses 
should not be regarded as small. Moreover, there has been considerable interest in kaonic 
excitations in dense nuclear matter leading, for example, to kaon condensation [^]- W% - 

Accordingly, it is useful to study strange excitations, and the strange quark mass is 
certainly not small on nuclear physics scales. The fact that the sum rules hold for any value 



of the quark mass was used explicitly in Ref. |12] in which the Nambu-Jona-Lasinio model 



was studied at mean-field level and in which the contribution to the sum rules coming from 
the zero-sound mode as well as from pion modes {i.e. modes that map smoothly onto the 
vacuum pion modes as the density goes to zero) were explicitly calculated. 

There is one major drawback to the sum rules derived in [jl],0: they only give information 
about modes with zero momentum. This is unfortunate because one of the most interesting 
issues about nuclear matter is the dependence of the energy of a mode on its momentum. 
The previous sum rules are completely insensitive to such effects. In the present paper, we 
will show how to exploit relativistic covariance to derive sum rules which are sensitive to 
states with nonzero momentum. More precisely, they involve sums over states with a fixed 
velocity, (3 = q/E, where g* and E are momentum and energy of the state, measured relative 
to the ground state of nuclear matter. 

We begin by stating the result in the rest frame of the medium; a covariant formulation 
will be given later. Also, for simplicity the notation is for two flavors. The generalization to 
three flavors is straightforward. Let us denote the ground state of translationally invariant 
nuclear medium by \C). The state is subject to space-independent constraints which fix 
its baryon density, pb, and isospin density, p/=i. We are concerned with the spectrum of 
excitations with quantum numbers of the pion (analogous results apply for the kaon) on top 
of the nuclear medium. These excited states are denoted by |ja,g), where ja labels modes 



with isospin equal (a = 0) , greater by one unit (a = +), or lower by one unit (a = — ) 
than the isospin of the state |C), and q is the momentum of the mode. The excitation 
energy of the mode, i.e., the difference of its energy and the reference energy of the state 
|C), is denoted by Ej^{q). We stress that no assumptions are made as to what the excited 
states \ja, q) are; they include collective excitations, one-particle-one-hole continuum, two- 
particle-two- hole continuum, etc. — in short, all states that have quantum numbers of the 
pion. 

We define the following spectral densities associated with the divergence of the axial 
current Af,: 

(^c{E,q) = 2^ ^1^ .-^1 5[E-EjSq)), a = 0,+,-. (1) 

In general, the sum is over continuum states and assumes the form of an integral. We will 

show that three classes of sum rules exist, which relate quark condensates and the isovector 

density to spectral integrals over a^: 

— f°° dE -> 

-mu{uu)c - md{dd)c = 2 -—a^{E,q = (3E) , (i) 
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_^ -^ (a+(E, q = m - ac{E, q = (3E)) , (iii) 
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Note that the spectral functions a^ in the integrands have the arguments constrained to 
q/E = (3. This means that the modes contributing to the sum rule move (in the rest frame 
of the medium) with the fixed velocity (3. 

Figure [^ illustrates schematically the described situation. What is shown is a typical 
result of a nuclear calculation of the charged pionic excitation spectrum in neutron matter 
up to the one-particle-one-hole level |T^. Solid lines correspond to the 7r+ and vr^ poles, and 
the shaded region represents the one-particle-one-hole continuum. In addition, a possible 
spin-isospin sound mode 7r+ [T^ is plotted. On the vertical axis is the frequency of the mode. 



uj. The energy of the excitation which positive (negative) isospin modes the energy is equal 
to + (— ) UJ. The dashed line shows the integration path in the sum rules. It corresponds to 



\(i\/E = (3, i.e., to modes moving with the fixed velocity (3 with respect to the medium (which 
is at rest). Contributions to the sum rules coming from various excitations are denoted by 
blobs (poles) and the thick line (the particle- hole cut). The dashed line is inclined to the 
vertical axis at the angle a = arctg/5. At various values of the velocity (3 different regions 
of the spectral density are sampled. Note that a < 45°. One should note that in reality 
the corresponding figure would be much more complicated. Due to multi-particle-multi-hole 
continua spanning the whole range of u and q, all states have finite widths, and contributions 
to sum rules are collected from everywhere along the dashed lines. 

We now pass to the proof of sum rules (|iD-(|i|). We use two facts. Firstly, the chiral 
current algebra satisfied by QCD yields the following operator identities: 

[Qt, [Qt, nm] = ^(0){r72, {r72, M}}g(0) , any a = 1, 2, 3 , (2) 

[g«,A^(0)]=ze"^%^(0), (3) 

where V^{x) = ^(a;)7^|r'^g(x) and A^(a;) = ^(a;)757^|r'^g(x) are vector and axial currents, 
Qs = Jd^xA^lx), 7i{x) is the QCD Hamiltonian density, and M = dmg{mu,md) is the 
quark mass matrix. 

Secondly, we assume that the medium is translationally invariant. This is true for infinite- 
volume nuclear matter. The important observation is, however, that the medium need 
not be at rest. Medium moving with a constant velocity —(3 in a reference frame is also 
translationally invariant. We denote such a state by \C, —f3). The sum rules are constructed 
in the usual way [jn[. The identities (0)-(0) are sandwiched by the state \C, ~j3). Inside the 
LHS we insert covariantly normalized intermediate states |I[], using the identity 

The intermediate states can be labeled by momentum p since \C, —(3) is translationally 
invariant. Index j sums over all additional quantum numbers. The excitation energy Ej (p) 
is the difference of the energy of the excited state and the state of the moving medium, i.e., 
we have 



-(/5)ir- _m H-u- ^ _ /p(/3)^ p(/3) 



H\C, -P) = E'a"\C, -P) , H\j,f^ = [E},'> + Ef\p)) |j,p} . (5) 

The momentum p is defined analogously — it is the momentum relative to the ground state of 
nuclear matter. Note that E- (p) and p form a Lorentz four- vector. The immediate result 
of the described construction is the following set of sum rules: 



{jo,p = 0\d-A\0)\C,-P)\[ 

0+,p = 0|9-A+(0)|C,-/3)|2 



io \Ej,{p = 0)\E)'^^'{p = 0) 



■{C,-/3\{mu + md)iuu + dd)\C,-P) = J2 .^(m,^ .Mr.(/3).- ^n 



\{j.,p = o\d-A-mc,-m' 

j- 









j+ \^j+ 



|(j_,p = 0|g-Ao-(0)|fi-/3)p 
J- \E)j[p = Q)V 

Sum rules ®-(0) are the consequence of Eq. (H), and sum rule (||) follows from Eq. @. 
The factor (1 — /3^)~^/^ is the dilatation factor for the isospin density, which is the time- 
component of a Lorentz four-vector (our notation is that p/=i is the isospin density in the 
rest frame on the medium). Note that only states with p = contribute to the above sum 
rules. 

We can now make a boost with velocity (3 to the rest frame of the nuclear matter. This 
boost transforms Eqs. (§[)-(§) into Eqs. (|i|)- ([!]). After this boost, the medium is at rest, 
and the excitations move with the fixed velocity 13. We can write the sum rules in the 
covariant form. Note there are two Lorentz vectors available for constructing invariants: the 
four-momentum p^ formed by the excitation energy and the momentum of the excited mode 
relative to the medium, and the four- velocity of the medium, u^. Introducing s = p^p^ and 
y = PiJ.u^ we can rewrite (|iD-(|iiI) as 

-mu{uu)c - mSd)c = ^H ^a^(s = y\l - (3^)\ y) , (I) 

Jo y 

-(m^ + m,){uu + dd)c = 2r ^(aMs = y\l-/3'),y) + acis = y\l-f3'),y)), (II) 

J —oo y 
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-{V^mc = uT % (a+(. = y\l - P%y) - a^s = y'il - P%y)) , (III) 



Covariant forms of the spectral densities (HD appear in the above equations. However, 
the arguments are constrained to s = y^ (1 — /?^) where /3 is the velocity of the modes in 
the medium's rest frame. An alternative — and equivalent — covariant formulation would be 
to use s as the integration variable in (|ID-(p^), with y constrained. Such a formulation is 
in some sense more natural in that it has an obvious vacuum limit — in the vacuum the 
spectral density is independent of y and the sum rules are written in terms of integrals 
over s. However, in the medium this is slightly awkward to do since negative energy states 
relative to the nuclear matter are possible — e.g.^ in neutron matter a state with vr''" quantum 
numbers may lie below the neutron matter ground state — and hence s can be a multi-valued 
function of E and the integral must be extended over each branch (see comment (4) below). 
A few remarks are in place: 
(1) If \C) is the vacuum, then (|)-(|i|) become the Gell-Mann-Oakes-Renner sum rule [jl5 



In the chiral limit and for vtiu = rrid = m one obtains the familiar relation Tn{qq) = F'^iv?^. 
Note, however, that the sum rules are exact for any values of uiu and 171^., also far away from 
the chiral limit. 

(2) Sum rule ( pi] ) or ( |III| ) is trivial for the vacuum, and also for isospin-symmetric nuclear 



matter. However, it is nontrivial for a medium which breaks the isospin symmetry. Its form 
is reminiscent of the sum rules of Fubini and Furlan (see e.g. [0]). 



(3) In the vacuum the sum rules for various values of (3 are equivalent. This is be- 
cause the dispersion relations for pionic excitations are fixed by Lorentz invariance, i.e. 



E = J q^ + m^. This is no longer true in the presence of the medium, and sum rules with 
different (3 are physically distinct. The point is that if one were to apply one of our finite 
velocity sum rules on a Lorentz-invariant state one would find that the size of the contribu- 
tions to the sum from any given mode would not depend on (3. In contrast, since the medium 
breaks Lorentz invariance one finds that the size of the contribution from a given mode to 
the sum does, in general, depend on the value of (3. In this sense these finite velocity sum 



rules provide additional information about the spectrum which is not present in the f3 = 
sum rules of Refs. [0]- 0. 

(4) The integration variable in the sum rule for the a = excitations, (|i|) or (j!]), ranges 
from to cxD, which reflects the fact that the state \C) is the ground state of the matter 
subject to constraints. Therefore all excitations within the constrained space have to raise 
the energy of the system. Hence, the integration variables in (|) and (|) are positive. This is 
not true for excitations with a = ±, which take the system out of the constrained space ^. 
In that case, states \j±, q) may have lower energy than \C), and the integration in (|ii|)- (H ) 
and (PI)- (PI) has to range from — oo to oo. A model example of such a behavior is given 
in ||12[. It is for this reason that it is awkward to write the covariant versions of these sum 
rules over the s variable. 

(5) In models where the pion fields satisfy the partially-conserved axial current condition, 
i.e. where d ■ A"" = —iv?^Ft^'k°' , the spectral densities (|I|) are proportional to the imaginary 
parts of the pion propagator. For such models the sum rules can be written as dispersion 
relations for the in-medium pion propagator. 

(6) The final remark concerns renormalization. Strictly speaking, sum rules (0)-(0) or (|)- 
(y) are ill defined on both sides. The left-hand side has an ill-defined composite operator and 
the right-hand side has a divergent sum. In order to make sense of the sum rules one needs 
to define some scheme to renormalize the qq operator, as well as a subtraction term for the 
spectral sums. At first blush this seems to suggest that these sum rules are useless. However, 
the need for renormalization stems from the vacuum sector of the theory. No new divergences 
are induced from the presence of the medium. Thus, once the renormalization is carried for 
the vacuum sector, it holds also for the medium. This means, that the sum rules (|)-(|ii|) 
and (|)-(||) should be regarded as vacuum-subtracted sum rules, involving the difference of 
the in-medium and the vacuum values of the quark condensates, {C\qq\C) — (vac|gg|vac), 
and accordingly subtracted spectral densities. The third sum rule, (|iIiD or ( |11I|) involves a 
conserved current, and as such requires no subtractions. 

To make things concrete, we will illustrate how our sum rules work out for a simple 
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toy model [|l[. The model describes the pion moving in an isospin-asymmetric medium and 
interacting with it only via p-meson exchange. Such a model is obviously quite unrealistic 
since, among other things, it assumes that the vr-N coupling constant, and hence qa is 
strictly zero. We set mu = irtd = m. In order to ensure chiral symmetry in such a model 
one must assume universal coupling of the rho meson and the KSFR relation |jl6|], 2g'l = 
m? / F"^. The inverse-charged pion propagator in the rest frame of the medium has the form 
G^{(l0i Q) = QoT Pi=i/F^Qo ~ (f ~ f^l- The model has the following spectrum of neutral, 
positive and negative isospin excitations: 



E+{q) 
E-{q) 






-pi=,l{2F^) + v/(p/=i/(2F^2))2 + ml + q^. 



(9) 



After some simple algebra, sum rules (|)-(|ii]) can be cast in the form 



-m{qq)c = F^ml 
—rn{qq) — ^^™^ 



)c 



^>; 



a-iP) 



+ 



a+iP) 



-Pi=i = -pi=i 



a_(/5) + a+(/?) a_(/5) + a+(/3) 



a_(/3)2-a+(/5)2 a^/?)^ - a+(/3)^ 



F^rnl , 



= -Pi=i 



(10) 
(11) 
(12) 



where a±(/?) = ±pi=i/{2F^) + ^ {pi=i/ {2F^)y + ml{l - P^). Sum rule (|) acquires the 
trivial form ([ToD . Positive (negative) isospin contributions to the sum rules (|i|) and (|iJD 
correspond to first (second) terms in the brackets in Eqs. (pJ] ) - ([T^) . The relative contribu- 
tion of the positive and negative isospin modes depends on the value of (3. These relative 
contributions to the sum rules are plotted in Fig. |^. The figure is drawn for negative pi=i 
(neutron matter). Note that as /3 changes, the relative weight of the modes changes. This 
illustrates our central point — for translationally-invariant systems which break Lorentz in- 
variance, these finite velocity sum rules are inequivalent to their zero-velocity counterparts. 
As /5 -^ 1, the positive isospin mode saturates both sum rules. Also note that the two modes 
saturate the sum rules (Eifl-fni), since there are no other excitation in this model. 
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Obviously, the toy model is not realistic, but it illustrates our general statements: distinct 
sum rules for different velocities (3 and the exactness of the sum rules for any value of 
m. A more realistic model would involve at least one-particle-one-hole-excitations. Cuts 
associated with such excitations contribute to sum rules in addition to poles. An example of 



such a model can be found in Ref. [0], where it was found that the contribution of cuts to 
the sum rules (at (3 = 0) was very small. In reality, as already mentioned in the discussion 
of Fig. m, there are in addition multi-particle-multi-hole continua; hence, all quasiparticles 
acquire finite widths. 

While the discussion heretofore has focused on pions in isospin asymmetric nuclear mat- 
ter, it should be obvious that the results go over mutis mutandus to the study of excitations 
with kaonic quantum numbers. The key point is that these sum rules work for any quark 
mass. Thus one can simply replace up quarks by strange quarks everywhere in the sum 
rules (or alternatively down quarks by strange quarks). The only change is that instead of 
the isospin density one will have the u-spin or v-spin density on the left-hand side of sum 
rule (iii) or (111). This is of significance since there has been considerable interest in kaons 
in nuclear matter and in the possibility of kaon condensation. 

It is worth noting that although these sum rules give information about the spectrum 
away from p = 0, they nevertheless do not make contact with all possible modes. The 
reason is that these modes are fixed velocity with (3 = p/E < 1 . Now one should notice the 
p/E is the velocity of the mode and not its group velocity. There is nothing in principle to 
prevent the existence of modes with j3 > 1. Our sum rules however tell us nothing about 
such modes. Recall, for instance, the p-wave pion condensation. If this were to happen, 
as one approached the transition density a mode with finite p would approach zero energy. 
Such modes do not contribute to our sum rules. 

Finally, we should discuss how our sum rules might be useful. Since they are not based 
on finite nuclei, it is hard to see what the sum rules can tell us about experiments directly. 
On the other hand, the sum rules do provide very strong constraints on model building. 
After all, these sum rules were derived from very basic properties known to be satisfied by 
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QCD — notably current algebra. Thus, any model which satisfies the various chiral Ward 
identities must satisfy our sum rules. Of course, it is very simple to construct models which 
satisfy the Ward identities. For example, mean-field models based on chiral langrangians 
will. However, such models typically exclude much essential nuclear phenomenology such as 
the effects of short-range correlations and the need to put in form factors to cut off spurious 
high-momentum physics. A formalism including such effects may well violate the Ward 
identities — that is to say, the approximations used to make the calculation tractable while 
including these effects will in general not be symmetry-conserving. One obvious use of our 
sum rules is to test how badly the symmetries are violated. 
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FIGURE CAPTIONS 

FIG. 1: Schematic plot of the excitation spectrum of neutron matter, involving 7r+ 
and TT^ poles, one-particle-one- hole continuum, and possible spin-isospin sound mode n^. 
Positive (negative) isospin modes have the excitation energy E equal to + (— ) uj. The 
dashed line shows the integration path in the sum rules, and contributions from various 
excitations are denoted by blobs (poles) and the thick line (cut). The dashed line is inclined 
to the vertical axis at the angle a = arctg/3. At various values of the velocity f3 different 
regions of the spectral density are sampled. Note, that a < 45°. In reality, due to multi- 
particle-multi-hole continua all states have finite widths, and contributions to sum rules are 
collected from everywhere along the dashed lines. 

FIG. 2: Toy model. Relative contribution to the sum rules as a function of velocity [3. 
The isovector density is set arbitrarily to p/=i = —F^rriT^. 
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